Exchange cotunneling through quantum dots with spin-orbit couphng 
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We investigate the effects of spin-orbit interaction (SOI) on the exchange cotunneUng through a 
spinful Coulomb blockaded quantum dot. In the case of zero magnetic field, Kondo effect is shown 
to take place via a Kramers doublet and the SOI will merely affect the Kondo temperature. In 
contrast, we find that the breaking of time-reversal symmetry in a finite field has a marked influence 
on the effective Anderson, and Kondo models for a single level. The nonlinear conductance can 
now be asymmetric in bias voltage and may depend strongly on direction of the magnetic field. A 
measurement of the angle dependence of finite-field cotunneling spectroscopy thus provides valuable 
information about orbital, and spin degrees of freedom and their mutual couphng. 
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Quantum dots based on materials with pronounced 
spin-orbit interaction (SOI), such as InAs, SiGe, carbon 
nanotubes, and single molecules have recently received 
reinforced attention.'^'Sl This is partially motivated by 
the quest for achieving electrical control of single spins, 
utilizing the fact that an electrical coupling to the or- 
bital degrees of freedom may allow for manipulations of 
the spin via the SOL^i^HlI! in quantum dots, the precise 
form of the spin-orbit coupling depends strongly on the 
band structure, confining potential and dot geometry al- 
together. It would therefore be of great value if one could 
infer about the SOI directly from a measured cotunnel- 
ing bias spectroscopy, which is known to produce sharp 
spectroscopic features due to threshold processes and/or 
Kondo effects. 

It is well known that Kondo effect in metals with mag- 
netic impurities like Ce and Yb, say, are strongly af- 
fected by spin-orbit interactiorflSHIll, The SOI modifies 
the spectrum of the impurity atom, but since it preserves 
time-reversal invariance a Kramers-degenerate ground- 
state remains and gives rise to Kondo effect. Likewise, 
a quantum dot holding a net spin-1/2, will also have its 
spectrum modified by SOI, and a Kramers degeneracy 
will still be available for Kondo effect. Unlike the atomic 
L • S coupling, however, the SOI in a quantum dot breaks 
rotational invariance and relates to specific spatial direc- 
tions, a kin to the effect of a crystal fields^*^"^ or nearby 
surface^i^'^ in the atomic problem. Since a spinful quan- 
tum dot allows for local directional probes such as bias 
voltage and magnetic field, the question arises if there 
are effects of SOI that show up directly in a transport 
measurement? 

Here we show that in the case where a single level 
approximation is valid, the SOI can be absorbed in a 
redefinition of the lead electron fields and thus leaves 
the Kondo effect unaffected. In the presence of a fi- 
nite magnetic field, however, spin, and orbital contents 
of the Kramers doublets become disentangled and a spa- 
tial asymmetry in the tunneling amplitudes can cause the 
Zeeman-split Kondo peak to become asymmetric in bias 
voltage. This type of asymmetric splitting does not occur 



without SOI, unless the voltage becomes large enough to 
allow for real charge fluctuations on the dot. Further- 
more, the SOI induced asymmetry will depend strongly 
on the direction of the magnetic field. The distinct an- 
gular dependence provides a very direct signature of the 
SOI in a quantum dot, thus providing valuable informa- 
tion about the quantum dot in question. 

We employ the following general single-particle Hamil- 
tonian to describe a quantum dot defined by a potential 
V{v) and placed in an external magnetic field 

Ud '- + V{v) + .g/isB • r 

+ -4^[E(r)x(p-eA)].r, (1) 

with r denoting the vector of Pauli matrices, and B the 
external magnetic field corresponding to a vector poten- 
tial A. The spin-orbit coupling is here kept on its most 
generic form in terms of the relevant nuclear or struc- 
tural electrical field E(r). The potential contains both 
the periodic potential from the ionic background and the 
imposed confining potentials defining the dot. 

In the absence of an external field (A = 0), this Hamil- 
tonian is symmetric under time reversal and its eigen- 
states therefore take the form of Kramers doublets of 
two spinors^-: 

where the wavefunction components u„ and w„ depend 
strongly on the confining potential. The corresponding 
eigenenergies, £„, come with a characteristic level spacing 
set by the confining potential and the strength of the SOI. 
Also the source, and drain electrodes may experience a 
SOI, so in general, we can express the eigenstates of the 
corresponding Hamiltonian in the leads, Hl/Rj in ths 
same way: 
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where a = L,R refers to left, and right lead, respectively. 
Using these eigenstates, the total many-body Hamilto- 
nian is given by 



■Dt = L/ R 

k. 



n,ri 



(4) 



a = L/R 

k, I/, ?7,n 

where c^^.^, creates an electron in the i^'th component of 
the Kramers doublet (z/ =t / J,) with momentum k in 
lead a — L/R, and dj^^ creates an electron in the ry'th 
component (77 / of the n'th Kramers doublet on 
the dot. For the interaction term we employ the con- 
stant interaction model Hint = Ec{N — Ng)'^, where Ec 
denotes the total capacitive charging energy of the dot. 

The amplitudes for tunneling between dot and leads 
t"^'" depend on the index of the Kramers doublets and it 
is given by the Hamiltonian overlap. 



= Jdr Ck.(r)Ht.t(r)V'„,(r) 



(5) 



where the total single-particle Hamiltonian still takes the 
form of Eq. ([T]) , but with an extended potential defining 
two tunneling barriers which support the distinction into 
leads and dot (jHtot = 'Rd + 'HL+'HR) made in our defini- 
tion of the eigenfunctions for the separate parts. Regard- 
less of the details of this potential, this first-quantized 
Hamiltonian takes the following form: 



Htotir) = Ho(r)rO + ^X,oe^JkT' Ej{r)d^ 



(6) 



with kinetic energy and local potential contained in 
7^0 (r) and the local spin-orbit term written in terms of 
the Levi-Cevita symbol Eijk (Einstein summation con- 
vention implied). Using the fact that the different 
Kramers doublet components can be related via time re- 
versaP3, i.e. 



(7) 



together with the relation r^T^ry = — (r')*, it is readily 
demonstrated that 



■'\<ii\Htot\n^) =-{a\<i^\Htot\ni^)\ 



(8) 
(9) 



which renders the tunneling ampHtude proportional to a 
unitary matrix in vq space: 



j.Qk:?T, J. TTjak? 



(10) 



fermion fields in the two leads: c 



Qk?7 



, TTjakm -p 



sufficiently large level spacings, we thus end up with the 
following single-orbital Anderson model: 



akr? ri 

~t~ ^ ^ ^gkm (^^akjj^^^V ^mr;^Q^kr;^ ^" Hint: (-^-^) 
ctkr; 

which no longer bears any trace of the SOI. Notice that 
this unitary transformation is specific to the m'th level 
and therefore tunneling amplitudes to any of the other 
levels on the dot will in general retain their full (unitary) 
2x2 matrix structure in vq space. Apart from its influ- 
ence on the precise magnitude of takm, SOI thus appears 
to have no effect whatsoever on transport phenomena in- 
volving only a single level. In particular, the Kramers 
degeneracy of this level will give rise to Kondo effect. 

This conclusion changes dramatically in the case of a 
finite applied magnetic field, which couples directly to 
the constituent quantum numbers of the Kramers dou- 
blets, i.e. to spin, and orbital degrees of freedom. Using 
symmetric gauge, A(r) = (B x r)/2, the magnetic field 
enters lid through the kinematic momentum. This gives 
rise to the following first quantized terms: 



'Hb=- MsB • L + MS ^(r^B^ - (r • B)^ 



(12) 



Mb 



3B + ((E(r).B)r-(E(r).r)B) 



Next, we consider a specific charge state with an odd 
number of electrons on the dot and assume all levels be- 
low the m'th level to be doubly occupied. For the singly 
occupied m'th level the dimensionless unitary matrix in 

Eq. (10) can be now be absorbed in a redefinition of the where electron creation operators, and cj^kn" ^^^^ 



with an orbital term depending on the angular momen- 
tum operator L, a diamagnetic term quadratic in i3, 
and a local anisotropic Zeeman term. The terms lin- 
ear in B both break the time-reversal symmetry and 
thus destroy the degeneracy of the Kramers doublets. 
We shall assume B to be weak enough that this split- 
ting, which we parameterize by an effective g factor, 
is much smaller than the relevant zero-field level spacing, 
i.e. q^ibB < Ae = min[e„+i - - £m-i]- 

Apart from this renormalization of the Zeeman split- 
ting within the m'th level, the linear terms in B also 
have off-diagonal terms which couple the state \mrj) to 
other states \nrj') via L and r. The amplitudes for tun- 
neling into the resulting finite B eigenstates of the dot 
are therefore changed and in particular the unitarity of 
^Qkm VlSqA iov B = Q is no longer guaranteed. In general, 
the matrix of tunneling amplitudes can be polar decom- 
posed into a product of a unitary and a hermitian matrix. 
The unitary part can again be absorbed in a canonical 
transformation of the conduction electrons in the corre- 
sponding lead and i^^m ^e taken to be hermitian 
in rj'rj space. Altogether, the tunneling term in (11) is 
therefore modified to 

Ht=Y: [ttJ.^.A^n + it'2j*dUc^^rr) , (13) 
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as tunneling amplitudes i^^„j now depend on the applied 
magnetic field. 

Kondo model: Within the Kondo regime, 

max(z^F,a|(C7l^)*i2,j"'|) < min(-£„,e™ + Ec), a 
SchriefFer- Wolff transformatiorp2i with the full r/'j] 
matrix tunneling amplitudes now leads to the following 
exchange-cotunncling (Kondo) model: 

= " Ma)5lkr,Cak^ + I^B9"^B'S^ (14) 



ak77 



9 '^a'a'-' ^a'k'77' 'i)'?7^ak)), 



ij—0,x,y.z 



with S'^ — 1, t"/^ — 6ri'n and cotunneling amplitudes: 



(15) 



where e± denotes the addition, and subtraction energies 
on the dot. Away from the particle-hole symmetric point, 
e+ = e_, a vector oi potential scattering amplitudes, J°-', 
is present. Notice also that expanding ( 15 1 to leading 

order in B, it follows from the hermiticity of i^^m ^^^^ 
the mira-lead exchange couplings will be diagonal and 
isotropic, i.e. f^,^ cx 

It is interesting to note that this exchange scattering 
has lead indices {L/R) mixed up with Kramers doublet 
indices (77) in such a way that the usual simplification 
to a single channel Kondo model no longer is possible. 
For 5 = (or without SOI) only one channel is in- 
volved, but a finite field breaks the L/R symmetry via 
the SOI and gives rise to a channel asymmetric two- 
channel (anisotropic) Kondo model. As we shall demon- 
strate below, certain system geometries will have zero 
Zeeman splitting in finite magnetic field and therefore a 
strong coupling two-channel regim^'-^'* should in fact be 
attainable for such geometries. 

Cotunneling current: From these cotunneling ampli- 
tudes one can now calculate the current through the dot 
apl 7 = eE„'^(r^,^ - r^«)P(77), The non-equihbrium 
occupation numbers P{ii) for the dot states satisfy a 
rate equation from which they are found to be ^(77) = 
^m/i^m + r^), with F^/,, = Y^a'a TJ^^'r?- These cotun- 
nehng rates are found as F^,^ = 7^, " Ae",^ns(Ae",,") 
where is the Bose function and the energy differ- 
ences are defined as Ae^,^ — im^ri' — £m,ri — Ma + 
fia' ■ Finally, the tunneling probabilities are 7"/," = 

real number since (jy,„ )* — jy„,. As for a system with- 
out SOI, the nonlinear conductance will exhibit cusped 
steps at bias voltage, V ~ ^ Ma, corresponding to 
the effective Zeeman splitting. Since, however, Kramers 
doublet and lead indices are mixed for finite magnetic 
field, the nonlinear conductance is no longer symmetric 
in bias voltage. In general, the two cusps at respectively 
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FIG. 1. Sketch of two-orbital model system characterized by 
spin-orbit field Bgo and the angular momentum vector A. 



positive and negative bias can now be of different magni- 
tude, and their relative magnitude will in general depend 
on the angle of the magnetic field. 

Two-level model: To better illustrate these results, we 
now exemplify our discussion by a simple two-level model 
(cf. Fig.[^. With two levels split by an energy S, we can 
express the Hamiltonian in the basis {|1 t); \^ i)i |2 t)i 
1 2 4,)}, where the wave functions of the two levels, ipi and 
"02, are chosen to be real. In this basis, the spin-orbit 
coupling is included to give the full dot Hamiltonian: 

{ia\Ho\ja') = -T^-rO,,<5/2 + rf^r^V Aso/2, (16) 

where we have chosen the spin quantization along the 
built-in spin-orbit field, fJ.s'Bso = 2^^(^il^('') ^ 
p\ip2) = — z i^so/"^, characteristic for these two levels. 
Hf) is diagonalized by two Kramers doublets: 



\ari) = u\lfi) 
\bri) = v\lri) 



IVT: \2Tfj, 



vv 



(17a) 
(17b) 



^so 



with u"^ + v"^ = 1, 2uv = Aso/A, A = ^/P 
and eigenenergies Ea/t — tA/2. Note that these dou 
blets follow the general structure of time-reversed pairs in 
Eq. ([2 ) . In the presence of a magnetic field, we shall ne- 
glect the quadratic term in ( 12 1 altogether. This amounts 



to assuming the dot to be much smaller than the mag- 
netic length, i.e. eBi^^f <C h, which will then ensure that 
eS2(r2)-e((r-B)2) 



^ ((rxp)-B). The remaining linear 
terms in Eq. ( 12 ) have the following matrix elements in 



the original four-state basis: 
{ia\HB\ja') 



goiJ-B ^0 p 

2 '-^ 



(18) 



where A = — i(l|L|2) is a (real) vector characteristic for 
the two levels. Expressing this in the zero-field eigenbasis 
(a, b doublets), we can now find the low field splitting of 
the individual doublets. That is, for fields low enough 
that this splitting will be much smaller than A, we obtain 
the effective single-doublet (m = a, b) Hamiltonians: 



{rari\H B\mri') 



9ol^B 



BZmism^rn,0,COS^rn)-Tnr,' (19) 



where we have defined 



Z. 



a/b 



5'^ sin^ C 



cosC ± 



2AA 



so 



90 A 



(20) 
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FIG. 2. (a) Geometry of vectors A, Bso, and B, indicating 
their relative angles. Bso and B span the x — z plane, (b- 
f) Nonlinear conductance, dl/dV, in arbitrary units vs. bias 
voltage in units of A for a set of representative parameters 
given in the insets. Each panel shows a progression of curves 
with varying ^, moving from the thick solid (black) curve 
with = to the thick dashed (black) curve with ( = n/2. 
Remaining parameters are qi-IbB — 0.5A, ip — 7r/4, ti,L/fl = 
1, t2,L = 3, t2,R = 0.1. Notice the closing of the Zeeman 
splitting for B _L Bso for two solely SOI-split orbitals in 
panel (c). This can also be seen analytically from Eq. (|20[). 



with sin^a/b — {sin ()/{Za/b^) and where A = 
|A|(cosCcos6' + sin ^ cos iy9 sin 0) denotes the projection 
of A on B in terms of the relative angles {C,9,(p) be- 



tween the two intrinsic vectors, B50 and A, and the ex- 
ternal B, all indicated in Fig. [2]ja). The eigenstates of 
this hamiltonian are readily found by a rotation within 
the plane spanned by B50 and B, i.e. linrj)^^'^ ~ 
Ryi^,n)\mr]) , where Ry{^) = cos(^/2) -I- ir^ sin(^/2). 
Apart from this splitting, a finite field will also mix 
the a and b doublets. Starting from this last eigenba- 
sis, which only refers to the direction of B, we include 
this mixing to linear order in the magnitude |B| from 
first order perturbation theory, i.e. lory)^^-' = lary)*^"^ + 
{bv'\HBH'myiUU"\bvT\ where {br,'\HB\av) = 
(XSt^,^ - go sin(C)AsoT^',,/2)MB-B/A. The amplitudes 
for tunneling into this lowest lying Zeeman-split Kramers 
doublet \ari)^^\ are now readily found using (17). Notice 
that it is only this last non-unitary mixing of a and b 
doublets which prevents us from diagonalizing t]!Jl by a 
canonical transformation of the conduction electrons. 

Whereas the SOI could not be discerned at zero field, 
the angle-dependent bias voltage asymmetry of dl/dV, 
confirmed by our simple model in Figs.[2jb-f), is a unique 
signature of SOI at finite field. Such bias asymmetries 
have ofte n be en found in experiments on various quan- 
tum dots.l^EH Nevertheless, it is often difficult to rule out 
the influence of incipient charge fluctuations, setting in at 
slightly higher bias voltages, as the source of this asym- 
metry. The only unambiguous evidence for such SOI- 
induced bias asymmetry will therefore be the observa- 
tion of its variation with a change in the direction of 
the magnetic field. Taken together with the possible an- 
gle dependence of the Zeeman splitting itself (cf. e.g. 
Refs. |3]and[5]), such a measurement can thus reveal oth- 
erwise inaccessible details on the spin-orbit coupling in a 
given quantum dot. 
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